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AVERAGE PROPERTIES OF COMPRESSIBLE LAMINAR BOUNDARY LAYER 
ON FLAT PLATE WITH UNSTEADY FLIGHT VELOCITY 1 


By Fkanklin K. Moore and Simon Ostrach 


SUMMARY 

The time-average characteristics of boundary layers over a 
flat plate in nearly quasi-steady flow are determined. The 
plate may be either insulated or isothermal. The time averages 
are found without specifying the plate velocity explicitly except 
that it is positive and has an average value. 

Each time average involves two groups of terms to the order 
considered in the report, a time average of quasi-steady terms, 
and terms related to the reduced frequency of the plate velocity 
fluctuations. 

The quasi-steady terms differ from the values for steady flow 
at the corresponding average velocity. These differences are 
reinforced by the frequency dependent averages for adiabatic 
wall temperature and heat-transfer rate. The effects oppose 
one another in the case of skin friction. 

The special case of harmonic velocity variation is considered, 
and it is found that large amplitudes accentuate the importance 
of the frequency-dependent terms. 

Oscillating the wall to increase the heat-transfer rale is not 
advantageous if the power to oscillate the plate is accounted for. 

INTRODUCTION 

In many current problems of aerodynamics, unsteady 
motions of a surface are important. The accelerating and 
decelerating phases of missile flight and the intermittent 
flow in an engine during unstable combustion are examples. 
The nature of the boundary layer of these unsteady flows 
may be studied with a view to determining friction drag, 
surface temperature, and rate of heat transfer through the 
surface. Usually, the boundary layer is so thin that' it 
responds almost instantly to temporal changes in flow con- 
ditions. Thus, the time history of the boundary layer is a 
succession of steady states, and such a boundary layer is 
called quasi-steady. If the motion involves accelerations 
which are particularly rapid, the quasi-steady description 
may require correction. 

In any case, it is often desired to assess the average effect 
of fluctuations in flow conditions. For example, if an insu- 
lated body is in motion with a speed that varies in time about 
some average value, the average friction drag and the aver- 
age surface temperature may differ from values appropriate 
to steady motion at the average speed. Another example is 
furnished by the speculation that the net performance of a 
heat exchanger could be altered by imparting an unsteady 
motion to the wall. 


In the present study, as an idealized special case of the 
foregoing type of problem, a semi- infini te flat plate is assumed 
to be in motion parallel to its surface and normal to its 
leading edge with a flight velocity U(t) that is always in 
the same direction but has a magnitude that fluctuates with 
time. The resulting boundary layer is assumed laminar and 
compressible, and the surface is either insulated or at con- 
stant temperature. The assumed Prandtl number is 0.72, 
which is appropriate for air under normal conditions. 

The boundary layer is assumed to be nearly quasi-steady. 
The basic boundary-layer analysis is already available for 
this problem in references 1 and 2, which treat the insulated 
plate and the constant plate temperature cases, respectively. 
The velocity profile in the boundary layer is found in the 
form 

g=| [F'(a)+UM+U'ffc)+ . . . -R UU<r)+ . • • ] 


( 1 ) 


where a is the usual Blasius variable signifying parabolic 
similarity in the boundary layer. (A complete definition 
will appear in a subsequent section.) The function F{a) 
is the Blasius function for steady motion of a flat plate. 
The parameters f n govern deviations from quasi-steadiness: 


X" & n U(T) 
i*- 1 £/»+! d2' n 


n — 1,2,3, . . . 


( 2 ) 


In cases of flight at substantia] speed, the factor U n+1 
usually ensures that f*_i is a rather small quantity. The 
functions F{<f), / 0 (<r), and fffo) are available in reference 1. 
A full list of symbol notation is provided in appendix A. 

If the plate surface is insulated, the temperature profile 
may be written as 

~-=l+ 2 ^ AU (2>(<r, f 0 , fi, h, • • •) (3a) 

where the “recovery factor” 

r=S(<r)4-f()ro( cr )4'firi(<r)+ . . . -HoTooMH- • • • (3b) 

and the functions R, r 0 , and r x are tabulated in reference 1. 
If the plate is at a constant temperature at which heat trans- 
fer takes place, reference 2 gives 


* Supersedes NAOA Technical Note 3888, "Average Properties of Compressible Tamlnar Boundary Layer on Flat Plato with Unsteady Flight Velocity,” by FranlUn K. Moore and 
Simon Ostrach, 18J6. 
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(4) 


and provides the functions H, S, ho, s 0 , hi, and Si. 

A study of the time-average properties of the flow repre- 
sented by equations (1) to (4), when the speed of the plate 2 
U fluctuates with time, is described herein. The functions 
/oo, r M , hoo, and are required for this purpose, and, not 
being available in references 1 and 2, are determined herein. 

A previous report (ref. 3) studies the average rate of heat 
transfer from an oscillatingjlat plate, hut differs from the 
present study in that the plate is doubly infinite and there 
is no net motion through the surrounding fluid; on the 
average, the plate is at rest. Appropriate comparisons are 
made herein between the. results of the present study and 
those of reference 3. 

DERIVATION OF FORMULAS FOR AVERAGE PROPERTIES OF 
BOUNDARY LAYERS 

The stream function p for nearly quasi-steady flows is 
defined in references 1 and 2 as 


where 


ip- — \]Cv a UX j(<r, f 0 , fi, f 3 , 

™ Uc 


.) 


u 


Ch m x 


(5) 


( 6 ) 


The coordinates X and T measured in a system with its 
origin fixed at the leading edge of the plate are related to ■ 
those (x,y) in a coordinate system which is stationary in 
the fluid by (see sketches (a) and (b)): 


X—x+^U df 

(7) 

Y=['f-dy 

(8) 


Equation (8) is emploj'ed to make the mpmentum equation 
independent of the energy equation. The velocity in the 
X,T-system is related to that in the 2 ,y-system by 


u(X, Y, T) =u a (x,y,t ) + U(t) 


(9) 


The constant O is the proportionality factor in the assumed 
viscosity-temperature variation 




e a 


=p 


C — 


( 10 ) 


’In the present problem, If compressibility Is Important, unsteadiness Is restricted to enter 
only through motion ol the plate. Thus, In a corresponding wind-tunnel test the model 
position would be varied mechanically and the tunnel flow would be held constant. The 
analyses of references 1 and 2 tor a compressible fluid would not apply for a tiled model In a 
pulsing flow. 


This constant may be evaluated by matching equation (10) 
with the Sutherland formula at some appropriate point, for 
example, at the wall. Thus 


„ IK /0„+216 o R\ 

V e a Wf 216° R ) 


(ID 


The unsteady boundary-layer characteristics can be deter- 
mined from equations (1) to (6). 

FORM OF UNSTEADY BOUNDARY-LAYER CHARACTERISTICS 

Insulated plate. — The unsteady boundary-layer charac- 
teristics from which the average properties will be determined 
are presented in this and the subsequent section in a more 
extended form than in references 1 and 2. The wall shear 

stress r„= (»a may be obtained in the following dimen- 
sionless form from equations (1), (8), and (10), and from the 
state equation for constant pressure (p0= Constant): 


° r ~ I 


5^ 


'gf [^"(0) + io/o(0) + f,/I(0) + ftToo(O) + . . .] 

( 12 ) 

The displacement thickness 5* may be defined 3 as 


(13) 


U(t) 


f U dt 


Position of ptole 
of lime /» O 


(a) 


um 


,-gexdY 


(b) 

(a) Coordinates fixed in fluid at rest, 
(b) Coordinates fixed in plate. 


> The steady-flow definition of displacement thickness is adopted herein. Actually, tbo 
displacement effect of the boundary layer is not properly represented by this definition If tbo 
flow is unsteady. However, In reference 4 It Is shown that the steady eipresslon Is part of tbo 
correct one, and because this expression yields the correct quasi-steady result Its calculation 
Is therefore warranted. 
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Upon application of equation (1), the state equation (j>9= Constant), and equations (3) into (13) there results 
d*=J2^^Slim(2a— F)— [f 0 /o(“)+fi/i(“)+lo/oo( o:> )+ • • [J 0 -^ d<r + fo J 0 r od<r+ 


fiJ o nd<r+fgj^ rood<r+ . . .JJ- 


where 

m=(y-l)M* (15) 

Isothermal plate. — The dimensionless wall-shear-stress, or local-skin-friction coefficients for the case of heat transfer at 
the plate, is identical in form with its counterpart for an insulated plate (eq. (12)). The constant C is evaluated by 
assuming that 0 W in equation (11) is the adiabatic wall temperature for the insulated plate and the maintained iso- 
thermal temperature in the heat-transfer case. 

The displacement thickness for this case is obtained from equation (1), the state equation, and equations (4) and (13) 
and is given by 

5*= { fan (2<r-F) - [fo/oO ) + fc/i( » ) +tf/oo( « ) + • • •]+ 


ml I iSdc-J-fo I Sodtr+fi 


where 


The local rate of heat transfer is given by . 


Q. 


r* a> r* cd —i 

v+fij Sidtr-H'oJ Soodtr-)- . . .J-(- 

24> [- Jo ^ ^odcr-|-f 1 J' hi do’-f-foj' At»d(r-j- . . .J^- (16) 

Therefore, for arbitrary fluctuations of the velocity, the aver- 
age local skin friction coefficient to order is obtained by 
applying equations (20), (22), (23), and (24) to equation (12) 
and integrating as in equation (21). Since for a periodic 
function, 

(18) , ! ^ ~° ( 25 ) 

Uol and also 

r r V==^- Cn+I) (^) 2 - (26) 


Using the definition of the Prandtl number, the state equa- 
tion, and equations (4), (8), and (10) in equation (18) yields following equation results: 


2— g: (0) s' <®+ 

to [>) 4-^ So(0) + . . -J+li [M(0)+ 


where 


^-hrrl • 

0 [m «o(0) + . . .]+r. [^(0)+ . |^(0)]p- s/ V) 2 } (27) 

— s ' ('o’) 4 - 1 I r' 1 r^(0) I m s' (01 I 1\ According to the remarks in the previous section, equation 

2$ 1 J 0 L 2$ 00 J J (27) holds for both the insulated and isothermal plates where 

(19) only the value of O differs. 

In a similar maimer the time averages of equations (14), 

VERAGES FOE ARBITRARY VELOCITY FLUCTUATIONS (16), and (19) 8X6 

assumed that the velocity fluctuates in a periodic m — — / 

he arbitrary manner, that is, /p^±. ( g -w u m (2a— F) f R dcr+ 

V U m \ <r->® Jo 

U(t) = U 1n g(cri) (20) rr , r «, 

the frequency of the fluctuations end „ is an J t nde+ 

isitive function so that />» % n 

, r3r roodAg-^igyl) (28) 

g=-Z-( g(r) dr=l (21) Jo y Jy 

0 for the displacement thickness of the insulated plate; 

T=wi (22) ¥=J^^(r^lim(2a-F)+m*^ f 0 Sda+ 

». t • /rt a \ 1 /a «-v\ • i . ... i • . /rt\ V %J 0 


TIME AVERAGES FOB ARBITRARY VELOCITY FLUCTUATIONS 

It is now assumed that the velocity fluctuates in a periodic 
but otherwise arbitrary manner, that is, 

U(i) = U K g(cd) (20) 

where w is the frequency of the fluctuations and g is an 
arbitrary positive function so that 


^/> )d - 


Substituting equations (20) and (22) into equation (2) 
yields 






with the frequency parameter O given by 

-Xco 


^ (g=& lira (2 a-F) J" S da+ 

dcr+Q 5 ^ 

8ld£r+ Jo ®«> dv^) g- bri (g'Y+ 

2$ + J\ da+ J o “ Aoo dc) g-^wry ) (29) 
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for the displacement thickness of the isothermal plate; and 

O s {[|«(0)+^(0)] g~™(gy+ 

<T 3 ' 2 G/') , \) ' (30) 

for the average local heat-transfer rate of an isothermal plate 
where 



(31) 


The average value of the adiabatic wall temperature (i. 
for the insulated plate) is also obtained in a similar manner 
from equation (3) and is given by 


SECOND-ORDER SOLUTIONS 

- In order to determine the average properties of unsteady 
boundary layers to order Q 2 , it is evident from equations 
(27) to (32) that the second-order (in f 0 ) solutions (denoted 
by double-zero subscript) must be .known. These solu- 
tions are determined in subsequent sections by extending the 
results of references 1 and 2. 


DIFFERENTIAL EQUATIONS AND BOUNDARY CONDITIONS 

The partial differential equations describing the unsteady 
flow and heat transfer of compressible viscous fluids are 
given in reference 1 as 

'f'rT'F'f'y'Pxr — 'i'x'pYY—U' (T)-\-0v< o \l/yyy (33) 

and 

Gv Cv 

Ot+tPyOx- ( try ) 2 (34a) 


For an isothermal plate, equations (37) are replaced (see 
ref. 2) by 


e(X,~,T)=e(0,Y,T)= 0 

(38a) 

Q(X,0,T)= 1 

(38b) 

SOLUTIONS 



The methods of solutions for the insulated and isothermal 
plate are identical to those in references 1 and 2, respectively, 
but, in each case, the method is extended herein to yield 
second-order results. 

Insulated plate. — For the insulated plate the boundary- 
value problem is defined by equations (33), (34a), (36), and 
(37). For nearly quasi-steady flows the stream function 
given by equation (5) and the temperature function given 
by equation (3a) are substituted into equations (33), (34a), 
(36), and (37) to yield 

U*+fU= -8^+2 ( 2 ^ +XpU M t; aX ) /*+ 2 ^ cf„- 

2X/ TO ^)/ r n tnX- S r&Yf'^]fct n f*r (39) 

n«=»0 U n-0 

r„+Pr/r,+y (j„Y=4Pr (2 ^H-§ 

< 40 > 

/»( co ,i'*)— 2; M0,f a )=/(0,f„)=0 (41) 

r,(0,r,)=r(®,f s )=0 (42) 

Introduction of equation (2) makes equations (41) and 
(42) self-consistent (i. e., functions of <r and only). For 
nearly quasi-steady flows, (see refs. 1 or 2) ; there- 

fore, the functions/ and r can be expanded as follows: 

/(<LfO =P(^)+fc/oM+fi/i.(<r) + . . . +fo/oo(v) + . . . + 
fbf/oi (<*) + • • • (43) 

=i2((r)-j-foro(<r)+firi((r)-l- . . . +fo r oo(v)+ ■ • • + 
fofi»oi ( 0 ) + ... (44) 


In reference 2 it was found more convenient to write equa- 
tion (34a) as 


where 


9 r+ 'PrQx— ifac© y — 0 yy+ ~~ (q Z. g y(&r 


r ) 2 


0 


e-e w 

@0 Qm 


(34b) 

(35) 


The appropriate boundary conditions on $ are 

Tpy(X,<x> ,T)=MO,Y ,T) —U(T) (36a) 

MX,O,T)=MX,O,T)=0 (36b) 

For the case of an insulated surface as treated in reference 
1, the boundary conditions on 0 are 


and 


e(X,o,T)=e(0,Y,T)=e a 


0y(X,O,T) = O 


(37a) 


Substitution of equations (42) and (44) into equations (39) 
to (42) and collection of terms independent of f and those 
multipled by f 0 and ft yield the three sets of equations solved 
in reference 1. Since for the present purposes the noxt 
higher order terms of equations (43) and (44) are needed 
(see eqs. (27) to (32), e. g.), the terms multiplied by fjj are 
collected and yield 

foo+F& — 4J?'f 00 -{-5F r foo = — 2/o(2 — _/o)-f/o(2cr — 3/ 0 ) (46) 

and 

^o+Pr(Fr' w — iF'r 0 o)=Pr[2or' a +2f ( ?o— SP'/oo— 

Hfo-Ffoo-^Y] (46) 

in addition to 

Joo( m ) = Joo(0) =/oo(0) = 0 (47) 

and 

r«(0)=r w (»)=0 (48) 

The function F and its derivatives are available in reference 


(37b) 
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5, and/o/o and their derivatives are presented in reference 1. 
Honce, the second-order functions and r M for the insulated 
plate are defined by equations (45) to (48). 

Isothermal plate. — Equation (45) holds for the isothermal 
plate as well as for the insulated one. However, for the iso- 
thermal plate, equations (34b) and (38) replace equations 
’ (34a) and (37) . The temperature function given by equation 
(4) is obtained by letting 


where 


e«=A0r,r«) 


U\T) 

2c p (e n —e a ) 


®( ff > ftx) 


and for the isothermal plate case, 

8*=(1.721)^^2 {(1+1.13 $)F I75 +(0.1673)mV^+ 

0* [(1.986+1.095 $) g- 9l *g'*+ (0.3112) mV 6/ y a J+ • • •} 

(56) 

2= (0.8211) (0.424) 

G J [(0.8350)§2- 7/ y i +(0.7574)mV 3 V il ]+ • • •} (57) 


^'(fl r )f#) =: if( <,r )+l'(A)(v)+fi^i(ff)+ • • • +fo*oo(«)+ • • • + 
?oftAoi(°") + • • • 

and 

s(<r,ft>) B <S r (tr)+foS 0 (<r)+fiSi(v)+ . . . +fo8oo(<t) + • • • + 
£ofi%( tr ) 

Substituting equation (4) into equations (34b) and (38) 
yields those equations treated in reference 2 independent of 
£ and to order £ 0 and ft. To order ft| the following equations 
are obtained: 

J&+Pr{Fh^-4F , h 0 o)=Pr(2vK-yX-5Hy oa + 

2hofo~~8hi)) (49) 

4+Pr(E4-4E'soo)=Pr ^2cai-SS'f M - 

2Uo+28X-F"fn-~ (/£*] (50) 

Aoo(“)=Aoo(0)=0 (51) 

s oo( OT ) = %i(0) — 0 (52) 

The fimctions on the right-hand sides of equations (49) and 
(50) are given in references 1, 2, and 5. 

Equations (45) to (52) thus define four boundary-value 
problems for the fimctions f w , r h^, and $oo* These prob- 
lems are solved for a Prandtl number of 0.72 by a numerical 
integration method described in references 1 and 2. The 
fimctions are presented in table I. 

RESULTS FOR AVERAGE PROPERTIES OF UNSTEADY 
BOUNDARY LAYER 

The following formulas are obtained by substituting the 
results of references 1 and 2 and those of the previous sec- 
tion into equations (27), (28), (29), (30), and (32), re- 
spectively: 

0,= (0-6640)^^^ GP- (1.306) fl 2 r^Y 3 + . . .] (53) 

For the insulated plate case, 

¥= (1.721) {F T75 + (0. 645) m*i^+ 

fl s [(1.986)r 9/ y s +(3.420)mV 6/ y a J+ • • •} (54) 

j^=l+ (0.8480) ^ [?+ (5.051) Q s ? r V s + . -J 


DISCUSSION OF QUASI-STEADY TEEM 

The leading terms (i. e., terms independent of Q) of equa- 
tions (53) to (57) reflect the nonlinear dependence of the 
physical quantities on U in a quasi-steady situation. Under 
the restriQtion that g is positive, all these leading terms are, 
of course, positive. 

Obviously; in the completely steady insulated case, when 
< 7 = 1 , all the leading terms in the braces equal 1. In the 
steady isothermal case, of course, only the factors involving 
g became unity. In the general quasi-steady case, when 
only the restriction that g=l is applied (eq. (21)), the mag- 
nitudes of the leading terms may be compared with unity 
(the steady case) using a special case of the Schwartz ine- 
quality (ref. 6), namely, 

' ~g**>zg m+n2 (58) 

If the choice m,«= 0,1 is made,- equation (58) yields 

0*251 (59a) 

For m, n=0, 1/2, the result is 

<+*:< 1 (59b) 

For m, 71=3/4, 1/4, g 312 g^>g= 1 and, therefore, in view of 
equation (59b), yields 

fpM (59c) 

If m=7i=l/4, <7 1/2 <7 _1/2 >1, which, together with equation 
(59b), yields 

^^>1 (59d) 

If m,7i=5/4,l/4, then p 73 g n *><F r ‘ a and, therefore, in view 
of equations (59b) and (59c), yields 

X 

g i,2 >l (59e) 

Insulated. surface. — Inspection of equations (53) to (55), 
together with equations (59), shows that_for_the insulated 
surface case, the quasi-steady values of O f , 5*, and 6 W are 
always greater than or equal to the corresponding quantity 
for unif orm flow at the average velocity U m . 

The foregoing effects are clarified by considering the aver- 
age velocity profile. For example, suppose that for half of 
the time, 0=1/2, and for the rest of the time, 2=3/2, so that 


(55) 
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<7=1 (fig. 1(a)); for incompressible flow, the average quasi- 
steady velocity profile is 



The result appears in figure 1 (b) ; the dash-dot lines represent 
the instantaneous quasi-steady profiles at the low and high 
velocities, the dash line represents the Blasius profile for the 
average (as if <7=1), and the solid line is the average profile 
from equation (60). The velocity from equation (60) is 
greater than the Blasius value near the wall and less than 
that value far from the wall. These profile differences are 
more pronounced if the variation of stream velocity is more 
extreme; if <7=1/2 for nine-tenths of the time and if <7=11/2 
for the rest of the time (fig. 1(c)), then 



The result appears in figure 1 (d) . The deviation of the aver- 
age profile from the Blasius shape clearly implies increased 
skin friction, and also suggests an increase of 3*. 

Note that the effects just discussed can be determined 
without fully specifying the function g. All that is needed 
is a specification of the proportion of time during which the 
various velocity values apply; in effect, a probability density 
distribution for velocity is sufficient. In deriving equations 
(60) and (61), the illustrative velocity functions were chosen 
for simplicity; of course, the present analysis would not 
apply near the sudden velocity changes which were postulated. 


Heat transfer. — If there is heat transfer at constant, wall 
temperature, the quasi-steady value of 3* from equation (66) 
is greater than the Blasius value unless $ is nearly at its 
maximum negative value, —1, which corresponds to a wall 
temperature of absolute zero. The heat-transfer rate q must 
be regarded as a function of $ and to* (eq. (66)) . The finer 
cross-hatched plane of figure 2 denotes g as a function of $ 
and to* for the average U„(g i=l). Positive q denotes heat flow 
out of the surface into the gas. If g is not always 1, then the 
first term of equation (57), independent of m*, tends to di- 
minish the magnitude of quasi-steady heat-transfer rate. The 
second term, proportional to m* and independent of <T>, 
always provides less heat flow out of the surface. The result 
is sketched as the coarser crosshatched plane of figure 2. 

Figure 2 indicates that variable stream velocity results in 
increased heat transfer to the wall only in the case of a cooled 
(2<0) wall, and even then only if to* has a value so that 


m*> - — — - — — 
-0.424 ^73_i 


(62) 


in cases of negative $. 
heat transfer is 


If $)>0, the condition of zero average 


* $ ff I/2 

0.424 gS/S 


(63) 


DISCUSSION OF FREQUENCY-DEPENDENT TERMS 

Equation (53) indicates that average skin friction is 
diminished in proportion to the squaro of the frequency. 
The dominant factor in this effect is the response of the 
boundary layer to rate of change of acceleration (/i(0) in 


. 

— Time ovf 
. Instanlon 
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1 
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:roge, u/U m 
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(b) 



o i o ^ i.o 

r/Z-t u/Um^u/U 

(a) Distribution of plate velooity with time. Moderate velocity 

variation. 

(b) Boundary-layer velocity profile. Moderate velocity variation* 

Figure 1. — Effect of plate velocity variation 
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, (c) 


(c) Distribution of plate velocity with time. Abrupt velocity variation. 

(d) Boundary-layer velooity profile. Abrupt velocity variation, 
profile of average velocity in boundary layer. 
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oq. (27)). With respect to skin friction, the quasi-steady 
term and the frequency term are in opposition. Subject to 
the requirement that Q«l, a flat plate mil experience 
a diminished average drag if 


r»2> ^-1 

(1.306) <r s/ Y 3 


(64) 


Owing to the appearance of g ' 1 in inequality (64), abrupt 
velocity changes would favor drag reduction. 

Insulated surface. — The average value of 5*, from equa- 
tion (64), will increase with SI*. The adiabatic wall tem- 
perature increases with Si 2 , which reinforces the effect of the 
quasi-steady term. 

Heat transfer. — From equation (57), the terms propor- 
tional to SI 2 all reinforce the effects of the quasi-steady terms 
with respect to the steady value. That is, the Q 2 term in $ 
reduced the magnitude of g as does the quasi-steady term 
since g ll *< 1, and the negative value of the Si 2 term in m* is 

Heal transfer 
| j Average 

| | Steady 



FiauRE 2. — Effect of unsteady velocity on relation among average 
heat transfer, temperature difference, and Mach number. 


con sistent with the corresponding quasi-steady term since 

^>1. 

In the heat-transfer case, the Q 2 terms of 3* also reinforce 
the corresponding quasi-steady terms, when the groups 
depending on $ and m* are considered separately, just as in 
the foregoing discussion of q_. 

RESULTS FOR HARMONIC ?(t) 

It may often be of interest to evaluate equations (63) to 
(57) when the function g is specialized to have harmonic 
form 

1+e sin r 

where the restriction 0<[5<T is imposed so that flight 
direction is not reversed. 

Arbitrary e. — The various integrals are evaluated for arbi- 
trary e in appendix B in the order of their appearance in 
equations (53) to (57) with the following results: 


[4E(*)-(l-e) #©] 


I%=g~ m g ,t = 


3vV 1+e 

2 K{k) 


[g-^(i)] 


ir-\/l+e 


(BIO) 

(B13) 

(Bl) 


/4 =; £r 8/ y 2 


lOSWl+^l-e 2 ) 2 

X[(3+16e-6e 2 ) #© -2(3+46-3^ K{k)\ (Bl5) 




V*=— 1- 


Vw 


(B4) 

(B5) 

(B2) 


Y (65) 


I 7 =g l!2 =^^E(k) 

7 T 

Is=7 ”=- [(23+9e2) EQc)-S{l-e)K{k)) (Bll) 


It— g 7/ V J = 


lit ) 53 tr 3/3 '<7 ,, = Vl+e 

7T 

where 


l±2f± 6e3 

lSirVl+efl— <?)L 1+e 


E{k)~^ 


E(k)—K(k) 


(B14) 

(B12) 


A=V2e/(l+.e) 


(B3) 


and K and E are the complete elliptic integrals of the first 
and second kind, respectively. The results are presented 
in table II. 

Small e. — In the event that e is very small, to a first 
approximation, each of the foregoing equations for arbitrary 
e may be replaced by making the appropriate substitution 
for N in the formulas 
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0^1+1 AT(lV-l) e 2 
g- V 2 =| 63 


(66) 


e Near unity. — If e is near unity, then k (eq. (B3)) is also 
near unity, or k'^l—k 2 is nearly zero; and the 

asymptotic relations 

K(k) «ln -pi — ; _E(&) « 1 

V 1—6 

apply (ref. 7). Equations (65) become 

^ 1 


in its own plane so that the relative stream velocity oscil- 
lates. The power required to oscillate the plate should 
properly he assessed to the system. Unless this power is 
less than the additional heat transferred, no heat is ab- 
stracted from the hot gas; also in this case that power could 
be used directly to increase the energy obtained with the 
surface fixed or in steady motion. Of course, if the oscilla- 
tions were inherent to the system, the assessment would not 
be required. 

If the power required to oscillate the plate is taken into 
account, the comparison is as follows: The excess of power 
required beyond that corresponding to steady flow is given 
by 4 


AP=£\m~-U n r u (U m )] dX=i dX 

~ (69) 




— - ( In 4 _ _ ^ 

9 9 1057T (1-eH yi 11- * 8 / 


-v 3. 


ST V 3 -* 


V2 


TIE .2V2, 


g 6 ^ 


2^ 

32V2 


<r 7/ V 3 -»- 


-v/2 1 


15s- 




3s- 

4 




V 1 — e 

(■“vfe" 2 ) J 


(67) 


Thus, as e—>l, equations (53) to (57) show that the effect 
of the Q, 2 terms is greatly accentuated because their coeffi- 
cients approach infinity from equations (67). For example, 
equation (53) becomes 


(7^(0.6640)^^^ [(1.20) -(0.3918) ~-+ . . .] (68) 

It should be noted that under these circ ums tances («— >1) 
the higher order terms in G 2 become significant. 


REMARKS ON HEAT-TRANSFER PROBLEM 

In a previous section it was shown that the heat-transfer 
rate with an oscillating surface differs from that with a 
surface in steady motion (see fig. 2). The relative merits 
of the oscillations, however, depend upon the particular 
configuration or application. For example, it has been 
determined that a greater heat-transfer rate to a cooled 
oscillating wall can be obtained for 4><0 provided that 
inequality (62) applies. This result suggests that, if one 
wishes to increase the rate of heat abstraction from a gas 
flow over a cooled surface (as, e. g., in a heat exchanger), 
it may be advantageous to oscillate the surface mechanically 


and the increment of heat transfer into the wall owing to the 
oscillation is 

AQ=~j*ft- q (U m ) ] dX (70) 

where X is the length of the plate. 

Equat ion (69) can be evaluated from equation (63), 
because gC f is given by the same formula as G f with the 
exponent of g increased by 1 wherever it appears in equation 
(63): 


AP= (0.6640) p.U* -^2 -j_^( 0 *)] (71) 

Equation (70) can be evaluated directly from equation (67), 
setting p = l to define q(U m ): 

AQ=(0.8211)c/„p« > U„^/^^[(1-F 5 )$+ 

(0.424)m*(p 73 — l)+^’(n J )] (72) 

The question is whether AQ^>AP. Therefore, equation 
(71) is subtracted from equation (72) to yield 

AQ-AP=(O.82n)c f 0 m pJJ^^^[(l-W 2 h- 

(0.385) (? 75 -l)m*+(<7(fl 3 )] (73) 

Thus, if AQ — AP)>0, 

^-l>(0.385)m*f^ (74) 

l—g m 


If inequality (74) is to hold, ~ — 1 must be positive. Now, 
the greatest value which 6jd a can approach consistent -with 

1 Tbe formula used for A P accounts for all the power required to maintain the plate In flight. 
If only that power required to oscillate the plate In a steady wind-tunnel flow or In a stationary 

heat exchanger were desired, the formula AP— (U— (IJr.ilX would bo used, which 

Is less than that given In equation (69) In the amount r.(l/»)] dX Equation 

(69) is selected for discussion on the grounds that the difference cited would usually bo obargod 
as an energy los3 of a practical airborne system. 
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heat flow into the surface is the insulated surface value from 
equation (55). Thus, inequality (74) becomes 

0.424mV>jp-l>0.385m*^= (75) 

or, 

— T 75 — 1 

It is doubtful that any reasonable function of g can be advised 
to satisfy inequality (76). Inspection of table II shows that 
inequality (76) cannot be satisfied by a harmonic function of 
any amplitude. 

For a stationary configuration, inequality (74) is computed 
using the A P expression in footnote 4 and again the inequality 
cannot be satisfied. Therefore, for air at normal conditions, 
the power required to oscillate the plate would exceed the 
extra heat transfer obtained by oscillation. 

The results of the foregoing discussion may be compared 
with the result of reference 3, wherein it was concluded that 
oscillation of a doubly infinite plate, in a fluid otherwise at 
rest would also result in an increased heat transfer to the plate 



CONCLUSIONS 

The unsteady laminar boundary layer on a flat plate in 
compressible flow has been analyzed for the case of time- 
variable velocity of flight with a view to describing the 
time-average characteristics of such a boundary layer. 
Flight velocity is assumed to vary slowly enough so that the 
resulting boundary-layer flow is nearly, but not quite, quasi- 
steady. The wall temperature has been assumed constant 
both along the plate and in time; further analysis would be 
required for the case of fluctuating surface temperature. 

In order to obtain time averages, the expansions of flow 

XU' (t) 

quantities must include terms in f 0 = — jji ’ fi 3 — ]j 3 > 

and fo- The terms for f 0 and ft were taken from previous 
work, while the fo terms were obtained by numerical integra- 
tion, and are presented in tabular form herein. 


The time-averages of skin friction and “displacement 
thickness” are presented, as well as heat-transfer rate at tfie 
surface, and for the special case of an adiabatic wall, the 
surface temperature. 

A significant amount of information may be obtained 
without specifying U(t ) beyond the requirements that U{i) 
remain positive and have an average value U m . Each time 
average involves two groups of terms to the order contem- 
plated in the present report: A time average of quasi-steady 
terms, and terms proportional to the inverse square of the 
characteristic time of the velocity fluctuation (i. e., the 
square of reduced frequency). 

The quasi-steady terms differ from the values for steady 
flow at the corresponding average velocity, owing to the 
nonlinear dependence of the physical quantities on U(t). 
In fact, especially for extreme variations of U(t) about the 
mean value, the average velocity profile is steeper near the 
wall and more gradual in its outer portion than the Blasius 
profile which applies at each instant. Thus, in the quasi- 
steady approximation, s kin friction, “displacement thick- 
ness,” and adiabatic wall temperatures are greater on the 
average than for the case of constant velocity. The magni- 
tude of the part of the heat-transfer rate that is independent 
of Mach number is less in the quasi-steady approximation, 
whereas the Mach number dependent part differs in the 
direction of less heat out of the surface into the gas. 

The differences cited are reinforced by the frequency 
dependent averages for adiabatic wall temperature and heat- 
transfer rate. The effects oppose one another in the case 
of skin friction. 

The various time averages are derived for the special case 
of harmonic velocity variation. Large amplitude affects 
chiefly the frequency-dependent terms, greatly accentuating 
their importance. 

The question discussed is whether it would be advantageous 
to oscillate the surface of a heat exchanger in order to take 
advantage of the increased rate of heat transfer to the wall, 
and it is concluded that the heat-transfer advantage would 
generally be vitiated by the power requirement for oscillating 
the surface against the action of skin friction. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 5, 1956 
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APPENDIX A 


SYMBOLS 


constant defined by eq. (11) 
local skin-friction coefficient 
specific heat at constant pressure 
related to stream function for flat plate in steady flow 
functions related to stream function for unsteady 
flat-plate flow, i=0, 1 , 2, 00, . . . 
function related to plate velocity, defined by eq. (20) 
temperature function related to steady isothermal 
flat-plate flow 

functions related to temperature for unsteady iso- 
thermal flat-plate flow, i=0, 1 , 2, 00, . . . 
thermal conductivity coefficient 
Mach number 

function related to Mach number, defined by eq. (15) 
constant defined by eq. (31) 
general exponent 
power 

Prandtl number 

total heat-transfer along plate 

local heat-transfer rate 

function related to temperature for steady insulated 
flat-plate flow 

functions related to temperature for unsteady in- 
sulated flat-plate flow, i=0, 1 , 2, 00, . . . 
function related to temperature for steady iso- 
. thermal flat-plate flow 

functions related to temperature for unsteady iso- 
thermal flat-plate flow, i=0, 1 , 2, 00, . . . 
time 

stream or plate velocity in Z-direction (see sketch 

(b)) 

mean velocity- 


'll. relative velocity in Z-direction 

u a absolute velocity in z-direction 

Z coordinate along surface measured from leading edge 

x coordinate along surface in system fixed in fluid 
(see sketch (a)) 

T coordinate defined by eq. (8) 

y coordinate normal to surface 

y ratio of specific heats 

5* displacement thickness 

e amplitude of velocity fluctuations 

dimensionless parameter, 7i=0, 1,2, ... .(eq. (2)) 

0 dimensionless temperature difference 

6 temperature 

ji absolute viscosity coefficient 

v kinematic viscosity coefficient 

p density 

tr dimensionless coordinate defined by eq. (6) 


$ 

* 

O 

0 } 


dimensionless coordinate 


u 


u„ 


Cv a x 

function related to time by eq. (22) 
local wall shear stress 
constant defined by eq. (17) 
stream function 

frequency parameter defined by eq. (24) 
frequency of velocity fluctuations 
Subscripts: 

w .evaluation at wall (F=0) 
oo , evaluation in stream (F— >°°) 

Subscript notation for partial differentiation is used when 
convenient. Primes denote ordinary differentiation. 
Superscripts: 

time average as defined in eq. (21) 


\ 
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APPENDIX B 

DETERMINATION OF INTEGRALS FOR HARMONIC OSCILLATIONS 


Given that g(r) = 1+6 sin r, the various averages appear- 
ing in the equations are to be found. The averages g~ in and 
g m can quickly be expressed as complete elliptic integais by 
replacing r by 2 z and taking account of the symmetry of 
the integrands: 

-7= T 75 S JL r d« _ 2 C Ifl d z _ 2 KQc) 

2ir J 0 A /l-fe sinr wVl + e Jo -y/l—Ji?sm 1 2 3 4 z sV l + e 

(Bl) 

pl=i- r* VH- 6 sin r dr=^i^ f'* jl-F sin* 2 dz 

W c/o 7 T Jo 


Thus, equation (B7) can be evaluated knowing equation 
(B6), which may be obtained from a recursion formula, as 
follows: From integrating by parts, 




4f>- T),,+S C0S ‘ T 4 t “7(2V+3>1"" (1+< 3 “ T) 

Xstar dT “iev+3) [/r <1+< ™ T)i,+idT “ 

J 2r (l+esinT)^dTj (B8) 


JbH±l E(Jc) (B2) 

TT 

where 

k^2e/(l+e) (B3) 

and KQc) and EQc) are complete elliptic integrals of the first 
and second kind, respectively (eqs. (773.1) and (773.3), and 
tables (1040) and (1041) of ref. 7). 

With help from equations (858.3), (436.00), and (436.03) 
of reference 7, these results are obtained: 

(l+esinr) 2 dr=l+^ (B4) 

(TV ,= lbJo ( 1 + esin T )~ 2 cos2t dr — ~“ 1 +7= : i CBS) 
The remaining averages are of the forms 

g N ^^lj\ 1 + e 8in T f*l dr (B6) 

and also 

g N+ *g /i J (l+e sin r) N+2 cos 2 r dr 

=— (1— €*)^J* (1+ 6 sin r) X+2 dr 

, 1 f 2r /1 I • si7+§ , 

-1 — I (l+eBinr) 2 dr 

^ Jo 

— ^ J (1 + e sin r) x+1 dr (B7) 


Combining equations (B7) and (B8) yields 

(^+|)(l-^/^-2(AT+2)/^+(iV+|y r+ 5=0 

(B9) 


Substituting equations (Bl) and (B2) into (B9) yields all 
necessary relations for evaluating equations (B6) and (B7): 


[4^-(i-6)X] 


^^Vl+e [(23+9^15-8(1-6)^ 

;) 

(£-*) 


g~ m g ' 2 


(B10) 

CBll) 

(B12) 

(B13) 


oTT-y/l + e 


o-o- ^+««-e^-2C3 + 4.-3<->X] 

(B15) 

The foregoing averages are presented .in table II as func- 
tions of e. 
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TABLE I.— SECOND-ORDER SOLUTIONS 



n 

/oo 

fco 

/» 

t 

rco 

rot 

/ 

h oo 

hot 

// 

«co 

tfco 

0 

1.08350 

0 

0 

0 

L 51298 

-a 13182 

0 

1.03328 

0 

.1 

1.05201 

.10728 

.00539 

-.45493 

L 48921 

13245 

-.01321 

.60292 

.08043 

.2 

.96756 

.20864 

.02126 

-.79778 

L 42573 

-.13329 

-.02649 

.31413 

.12530 

.3 

.84501 

.29952 

.04677 

-1.04778 

L 33275 

-.15283 

-.03981 

.12345 

.14660 

.4 

.69821 

.37683 

.08071 

-L 22022 

1.21876 

-.12966 

-.05296 

-.00100 

.16217 

.5 

.53920 

.43876 

.12102 

-L 32738 

L 09088 

-.12273 

-.06562 

-.08207 

.14773 

.6 

.37789 

.48460 

.16792 

-L 37922 

.95513 

11146 

-.07737 

-. 13533 

.13668 

.7 

.22196 

.51452 

.21801 

-L 38395 

.81660 

-.09577 

-.08776 

-.17062 

.12125 

.8 

.07702 

.52936 

.27032 

-L 34855 

.67967 

-.07608 

-.09638 

-.19353 

.10297 

.9 

-.05306 

.53042 

.32342 

-L 27929 

.54802 

-.05328 

-.10286 

-.20673 

.08287 

L 0 

—.16583 

.51932 

.37600 

- L 18215 

.42474 

-.02857 

-.10696 

-.21111 

.06160 

LI 

-.25980 

.49788 

.42694 

-L 06312 

.31232 

-.00329 

-.10855 

-.20687 

.04093 

L 2 

-.33452 

.46800 

.47530 

—.92844 

.21264 

.02116 

-.10764 

-. 19421 

.02081 

1.3 

-.38974 

.43162 

.62032 

—.78455 

.12694 

.04352 

10439 

-.17382 

.00237 

1.4 

-.42593 

.39069 

.56147 

-.63797 

.05582 

.06272 

-.09905 

—.14709 

-.01373 

l.a 

-.44399 

.34704 

.59837 

-.49498 

-.00077 

.07802 

-.09199 

-. 11611 

-.02692 

1.6 

-.44530 

.30245 

.63085 

-.36120 

-.04348 

.08895 

-.08360 

-.08334 

— . 036 S 9 

1.7 

-. 43189 

.25847 

. 65 S 88 

—.24104 

-.07347 

.09541 

-.07435 

-.05113 

-.04361 

1.8 

-.40028 

. 21648 

. 6 S 261 

-.13706 

—.09225 

.0753 

-.06466 

-.02166 

-.04721 

1.9 

-.37091 

.17753 

.70228 

-.05322 

-.10164 

.09577 

-.05496 

.00333 

-.04809 

2.0 

-.32905 

.14254 

.71825 

.01202 

—.10354 

.09073 

-.04562 

.02295 

-.04674 

2.1 

-.28380 

.11182 

.73092 

.05906 

—.09934 

.08314 

-.03600 

.03695 

-.04369 

2.2 

-.23782 

.08580 

.74078 

.08080 

-.09227 

.07381 

-.02905 

.04565 

-.03953 

2.3 

-.19386 

.06419 

.74823 

.10679 

-.08234 

.06352 

-.02216 

.04973 

-.03472 

2.4 

—.15345 

.04688 

.75374 

.11281 

—.07128 

.05298 

-.01635 

.05009 

-.02969 

2.5 

-.11804 

.03339 

.75776 

.11005 

-.06005 

.04278 

-.01166 

.04770 

-.02470 

2. 6 

-.08819 

.02304 

.76051 

.10287 

-.04934 

.03338 

-.00775 

.04347 

-.02022 

2.7 

-.06387 

.01559 

.76247 

.09172 

-.03959 

.02507 

-.00486 

-.03828 

-.01613 

2.8 

-.04494 

.01008 

.76369 

.07900 

-.03105 

.01799 

-.00271 

.03263 

-.01258 

2.9 

—.03052 

.00642 

.76452 

.06587 

-.02380 

.01225 

-.00121 

.02698 

-.00961 

3.0 

-.02006 

.00392 

.76507 

.05335 

-.01785 

.00776 

-.00022 

.02171 

-.00717 

3.1 

-.01271 

.00223 

.76528 

.04207 

-.01309 

.00440 

.00038 

.01704 

-.00523 

3.2 

-.00773 

.00134 

.76557 

.03233 

-.00939 

.00200 

.00069 

.01304 

-.00374 

3.3 

-.00448 

.00060 

.70555 

.02424 

-.00057 

.00040 

.00081 

.00976 

-.00260 

3.4 

— .0023 S 

.00038 

.76566 

.01774 

-.00449 

-.00050 

.00080 

.00714 

-.00177 

3.5 

-.00110 

.00018 

.76573 

.01268 

-.00297 

-.00111 

.00071 

.00509 

-.00116 

3.6 

-.00046 

.00006 

.70562 

.00886 

-.00191 

-.00133 

.00069 

.00355 

-.00073 

3.7 

-.00015 

.00013 

.76579 

.00604 

-.00117 

-.00133 

.00046 

.00243 

-.00044 

3.8 

.00009 

-.00003 

.76561 

.00403 

-.00067 

-.00124 

.00033 

.00164 

-.00022 

3.8 

.00010 

.00014 

.76575 

.00263 

-.00034 

-.00106 

.00022 

.00108 

-.00009 

■(.0 

.00019 

.00006 

.76576 

.00169 

-.00013 

-.00087 

.00012 

.00070 

-.00001 

4.1 

.00021 

.00013 

.76565 

.00106 

.00000 

-.00070 

.00003 

.00045 

.00004 



r to 


CO 


r* oo 


CO 




I /oo d<jT=>2.5490 

roo dir -1.0422 

1 A 00 tier 

0.1033 

roo d<r-0.0713 



Jo 


JO 


Jo 


Jo 



TABLE IL— EVALUATION OF INTEGRALS 


€ 

h 

it 

it 


h 

u 


/. 

u 

/io 

0.2 

1.00748 

a 02087 

1.00769 

0.02269 

1.020 

0.02059 

a 99747 

1.03747 

a 02168 

0.02050 

.4 

1.03019 

.09673 

L 03299 

.13704 

1.080 

.09107 

.98960 

1.14963 

.11282 

.08680 

.6 

1.06 S 72 

.29227 

1.08539 

.71028 

1. 180 

.24999 

.97522 

1.33554 

.43430 

.22032 

.8 

1. 12432 

.96153 

L 20016 

7.28120 

1.320 

.66663 

.95129 

1.69356 

2,41388 

.49770 

.9 

1.15938 

2.38511 

•1.33236 

63.68454 

L 405 

1.29413 

.93287 

1.74880 

10.88546 

.79898 

.95 

1.17904 

5.30960 

L 47463 

629.649 

1. 45125 

2.20252 

.92025 

1.83273 

46.84304 

1.10872 





